Application of Commercial Non-Dispersive Infrared 
Spectroscopy Sensors for Sub-ambient Carbon Dioxide 

Detection 

Michael J. Swickrath* and Molly Andersoir 

NASA Johnson Space Center, Houston, TX, 77058 

Summer McMillirr 

Engineering and Science Contract Group - Jacobs Technology, Houston, TX, 77058 

Craig Broermard 

Engineering and Science Contract Group - Hamilton Sundstrand, Houston, TX, 77058 


Monitoring carbon dioxide (CO 2 ) concentration within a spacecraft or spacesuit is crit- 
ically important to ensuring the safety of the crew. Carbon dioxide uniquely absorbs light 
at wavelengths of 3.95 micrometers (/im) and 4.26 /im. As a result, non-dispersive infrared 
(NDIR) spectroscopy can be employed as a reliable and inexpensive method for the quan- 
tification of CO 2 within the atmosphere. A multitude of commercial off-the-shelf (COTS) 
NDIR sensors exist for CO 2 quantification. The COTS sensors provide reasonable accuracy 
as long as the measurements are attained under conditions close to the calibration con- 
ditions of the sensor (typically 21.1°C (70.0°F) and 1 atmosphere). However, as pressure 
deviates from atmospheric to the pressures associated with a spacecraft (8.0—10.2 pounds 
per square inch absolute (psia)) or spacesuit (4.1— 8.0 psia), the error in the measurement 
grows increasingly large. In addition to pressure and temperature dependencies, the in- 
frared transmissivity through a volume of gas also depends on the composition of the gas. 
As the composition is not known a priori, accurate sub-ambient detection must rely on 
iterative sensor compensation techniques. 

This manuscript describes the development of recursive compensation algorithms for 
sub-ambient detection of CO 2 with COTS NDIR sensors. In addition, the source of the 
exponential loss in accuracy is developed theoretically. The basis of the loss can be ex- 
plained through thermal, Doppler, and Lorentz broadening effects that arise as a result 
of the temperature, pressure, and composition of the gas mixture under analysis. This 
manuscript provides an approach to employing COTS sensors at sub-ambient conditions 
and may also lend insight into designing future NDIR sensors for aerospace application. 


‘Analyst, Crew and Thermal Systems Division, 2101 NASA Parkway, EC211, Houston, TX, 77058, AIAA Member, 
t Analysis Lead, Crew and Thermal Systems Division, 2101 NASA Parkway, EC211, Houston, TX, 77058, AIAA Member, 
t Project Engineer, EVA and Health Systems Group, 2224 Bay Area Blvd., Houston, TX, Member AIAA. 

§ Project Engineer, CxP and Advanced Systems Group, 2224 Bay Area Blvd., Houston, TX, Member AIAA. 


1 of 14 


American Institute of Aeronautics and Astronautics 



Nomenclature 


co 2 

= Carbon dioxide 

COTS 

= Commercial off-the-shelf 

CRDS 

= Cavity ring-down spectroscopy 

ECLSS 

= Environmental control and life support 

EVA 

= Extra-vehicular activity 

h 2 o 

= Water 

IR 

= infrared 

ISS 

= International Space Station 

n 2 

= Nitrogen 

NDIR 

= Non-dispersive infrared 

0 2 

= Oxygen 

ppm 

= Parts per million 

psia 

= Pounds per square inch 

VSA 

= Vacuum swing adsorption 


I. Introduction 

C ARBON dioxide produced through respiration can accumulate rapidly within closed spaces. If not man- 
aged, the crew’s respiratory rate increases, headaches and hyperventilation occur, vision and hearing are 
affected, and cognitive abilities decrease. Consequently, development continues on a number of CO 2 removal 
technologies for human spacecraft and spacesuits. Terrestrially, technology development requires precise 
performance characterization to qualify promising air revitalization equipment. On-orbit, instrumentation 
is required to identify and eliminate unsafe conditions. This necessitates accurate in situ CO 2 detection. 

A variety of techniques exist to quantify CO 2 in the atmosphere. One low-cost yet sensitive technique 
is non-dispersive infrared (NDIR) spectroscopy. Infrared absorption is associated with vibrational and ro- 
tational energy transitions of covalent bonds that are susceptible to energetic perturbations at a particular 
wavelength. 1-4 A bond stretching mode exists at 4.26 micrometers (/mi) for CO 2 . Excitation leads to rota- 
tional energy transitions producing a multitude of distinct absorption lines. As a result, the Beer-Lambert 
law can be used to correlate an intensity change to gas concentration, A' = 1 — exp (— k v L ), where A' is the 
absorbed fraction of the light ( A ' = A/A Q where A 0 represents gas absorbance when the sample is devoid 
of the species of interest), L is the optical path length, and k v is the absorption coefficient. A number of 
commercial off-the-shelf (COTS) NDIR sensors are readily available for CO 2 detection. 

However, the accuracy of COTS sensors are highly application specific. Absorption lines undergo broad- 
ening through a variety of mechanisms that diminish the sensitivity and precision of these sensors. 5 First, 
spectral lines have a contribution to width which arises from the Heisenberg uncertainty principle as it applies 
to spectroscopic measurements. Specifically, excited states have a finite lifetime governed by a decay process 
such that absorption measurements produce ensemble-averaged results. Second, random collisions perturb 
the energy levels of the components of a gas sample. This is referred to as ‘pressure-broadening’ and is 
strongly influenced by system pressure where high pressure increases the likelihood of molecular interaction. 
Pressure broadening is also composition dependent since kinetic excitation is influenced by the momenta of 
the molecules involved in an excitation event. Lastly, thermal energy imparts random Doppler velocities 
to each molecule which serves to ‘blur’ a line over a spectral range. This indicates that detection is also 
influenced by system temperature. 

A number of pressures, under consideration for space exploration, are summarized in detail elsewhere. 6 
The International Space Station (ISS) and Space Shuttle were designed for 14.7 pounds per square inch (psia), 
and the Shuttle can also operate at a pressure of 10.2 psia. For extra- vehicular activity (EVA), pressures 
of 4.3 psia are under consideration while the spacesuit must operate at 8.0 psia to interface with vehicles 
and as high as 23.1 psia under contingency operations such as decompression sickness mode. Moreover, 
atmospheric composition is dictated by the Environmental Control and Life Support System (ECLSS) and 
can range from low to high relative humidity with several orders of magnitude higher CO 2 partial pressure 
than exist on Earth. Technology development necessitates hardware testing over a wide variety of pressures, 
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temperatures, and compositions. COTS-NDIR sensors are attractive from a testing stand-point as they are 
low cost and resolute; however, detection over a wide range of pressures, temperatures, and gas compositions 
poses a challenge to employing these sensors. This manuscript establishes one approach to the tempera- 
ture/pressure/composition compensation of COTS-NDIR sensors. Furthermore, theoretical techniques to 
generate mixed gas spectra as a function of temperature and pressure are developed for potential application 
in aerospace sensor design. 


II. Scaling Heuristics 

COTS-NDIR sensors provide an inexpensive and flexible means to quantify CO 2 concentration in a gas 
sample. These sensors rely on the approximation that absorption measurements under one set of conditions 
(i.e. test conditions) can be corrected through scaling based on another set of conditions ( i.e . calibration 
conditions). This method is referred to as the non-overlapping line approximation 7 and it tends to produce 
reasonable results as long as the pressure and temperature of a system does not deviate significantly from 
calibration conditions. 


A 

P 



(1) 


In Eq. (1), A represents the total absolute absorption, P is the total pressure, and u is the absorber 
amount (u = piL) which depends on the component density of the gas sample, pi [pi = yiP/RT assuming 
ideality where y 2 is the component mole fraction). The parameter L is the path length of the optics of the 
sensor. Consequently, the molar fraction of a component can be quantified via Eq. (2) through inversion of 
function /. 
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Using the non-overlapping line approximation, scaling is accomplished presuming yipi = t/, )0 Pi,o so that 
compensation is accomplished according to the following relationship where the function g is the inverse of 
/• 


T 

Vi = yi,o7fT9 
-L o 


(3) 


In Eq. (3), t/j >0 is the mole fraction reported by the sensor and P a and T 0 are the sensor calibration pressure 
and temperature. The mole fraction is proportional to absorbance which is also proportional to the voltage 
output of a sensor. These proportionalities are all implicitly captured within the coefficients of function 
g. Literature exists suggesting g can take the form of a polynomial. 1,2 Alternatively, in this work, g has 
been represented as a power-law expression. This approach is physically more representative of the physics 
captured by the Beer-Lambert law. It will be shown within this manuscript that the power law produces 
reasonable agreement with experimental data. 

As previously mentioned, not all gases have the same capacity to cause pressure broadening. This renders 
NDIR detection composition dependent. The composition dependency is generally accounted for through 
the notion of ‘equivalent pressure’, P e . The equivalent pressure modifies the reported value for composition 
specific spectral broadening effects. 1,2 Coefficients (f>i represent the broadening of spectral lines relative to 
nitrogen (i.e. <p N2 = 1). 


Pe = P 


Components 

1 + - !) Vj 
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(4) 


P e is substituted for P in Eq. (3). The typical major constituents of the gas-phase in space exploration include 
CO 2 , H 2 O, N 2 , and O 2 . However, the sensor compensation routines developed within this investigation 
were established based on test data in an oxygen-free environment. As a result, Eq. (5) represents the 
compensation function established in this work for carbon dioxide detection with COTS-NDIR sensors. 
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In Eq. (5), a and (3 are power-law coefficients accounting for pressure-broadening, </> H2 o is the foreign-gas 
broadening coefficient associated with water-vapor, and <fi C o 2 is the self-broadening coefficient for CO 2 . 

Expression (5) provides a closed-form analytical expression for off-calibration compensation of COTS- 
NDIR sensors providing the four compensation parameters can be determined. However, the complication 
arises that the actual carbon dioxide concentration depends on a function directly dependent on the carbon 
dioxide concentration. Consequently, expression (5) must be solved iteratively. This was achieved through 
a Newton-Raphson scheme in which an objective function, F, was defined along with a corresponding 
derivative with respect to carbon dioxide fraction F' = dF / dy C02 . 
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The Newton-Raphson method is a relatively efficient recursive technique seeking roots of a real -valued 
function with a property of quadratic convergence. Employing expressions ( 6 ) and (7), the Newton-Raphson 
scheme proceeds from iteration n to n + 1 as follows. 


Veil = Vco 2 


F 

F' 


( 8 ) 


After each iteration, a convergence check is performed r = — Uco 2 \ ■ When the successive error between 

iterations r, is below a pre-defined threshold, convergence is presumed. The actual fraction of CO 2 for the 
first iteration was set to the reported fraction of CO 2 per the non-compensated sensor reading to initialize 
the routine, = Vco 2 ,o ■ The tolerance was set to 1 x ICC 8 and convergence was generally achieved within 

three to four iterations. Situations where convergence required significantly more than three to four iterations 
did arise and were indicative of instrumentation errors ( e.g . sensor-drift requiring a re-zero and re-span 
procedure). Consequently, counting iteration during the recursive scheme provided an additional measure 
to identify and eliminate sensor anomalies real-time during testing. The iterative technique described above 
provided the methodology to correct COTS-NDIR sensors across a wide variety of temperatures, pressures, 
and compositions. As a result, this enabled the application of these sensors across a wide variety experiments 
performed at the NASA Johnson Space Center (JSC). 


III. Experimental Characterization 

The first step in verifying the iterative technique involved building a custom sensor test stand. A mul- 
titude of COTS-NDIR sensors were evaluated during this investigation from two primary vendors. Vaisala 
HUMICAP® HMT-334 relative humidity sensors and CARBOCAP® GMT-221/222 CO 2 sensors (Vaisala, 
Vantaa, Finland) were employed for the initial investigation. Subsequent investigations relied on the Li-Cor 
840A Gas Analyzers (Li-Cor Environmental, Lincoln, NE). Sensors were ranged from 0% to 1% up to 0% to 
5% C0 2 with respect to standard temperature and pressure. 

A custom test apparatus was built for the preliminary investigation to characterize sensor performance 
at a reduced pressure environment. The initial test stand was capable of regulating pressure and introduc- 
ing precise quantities of a dry pre-mixed CO 2 /N 2 gas source. A chamber with laboratory air was initially 
de-pressurized to near vacuum and then refilled with dry N 2 for a number of times to remove all poten- 
tially NDIR sensitive gas species. After the last evacuation procedure, the chamber was then refilled with 
the pre-mixed gas source to a desired pressure ranging from 2-15 psia. Once the pressure was achieved, 
the mass flow controllers were powered off and sensor measurements were collected for several minutes. 
Following this procedure, the mass flow controllers were actuated again until the next desired pressure 
was reached. This process was followed up until ambient pressure was achieved. A Rosemount Analyti- 
cal XStream HASXEE-IM-HS gas analyzer (Emerson Process Management, Pittsburgh, PA) was used to 
provide actual gas concentrations. 

These preliminary experiments provided a basis for the calibration function /. However, without the 
capability to introduce precise amounts of water vapor, complete compensation could not be achieved on 
this test-stand alone. 


4 of 14 


American Institute of Aeronautics and Astronautics 



Subsequent investigations relied on a pre-existing test stand established to collect performance data for 
amine swing beds being tested at the JSC. 8,9 This test stand was capable of precise CO 2 and water vapor 
introduction through a Bronkhorst Controlled-evaporator and mixing unit (Bronkhorst, Bethlehem, PA). 
The loop was de-pressurized for sub-ambient testing as low as 3.8 psia (pounds per square inch, absolute) 
using a Varian TriScroll™ 300 vacuum pump (Varian, Inc., Agilent Technologies, Santa Clara, CA). Gas 
flow within the loop was controlled using a Micronel® U51DX-024KK-5 fan (Micronel® AG, Tagelswangen, 
Switzerland) which was subsequently measured by a Teledyne Hastings HFM-200 flow meter. The test article 
allowed for the collection of transient data over a range of relative humidity and CO 2 partial pressures. A 
cavity ring-down mass spectroscopy (CRDS) unit (Picarro, Sunnyvale, CA) was plumbed into the loop in 
series with the CO 2 sensors. CRDS takes multiple light intensity readings as it travels between opposing 
optical mirrors over the course of micro-seconds. The intensity decays as the light traverses the optical 
path over several measurements providing data with which to establish precise species concentrations. The 
cavity within the CRDS is maintained at a constant pressure of 140 millimeters of mercuary (mm Hg) and 
10 °C (50 °F) irrespective of the pressure and temperature of the test stand. Consequently, the CRDS 
unit was considered as the actual reading providing data with which to find the appropriate constants in 
expression (5). Whereas the CRDS collects measurements over a wide temperature/pressure environment in 
which we place high confidence, the cost and size of the CRDS unit relative to COTS-NDIR sensors renders 
it prohibitive for wide-scale application. 

IV. Verification of the Compensation Method 

As previously mentioned, a custom test stand was built to explore the relationship between pressure and 
carbon dioxide concentration on reported readings for the COTS-NDIR sensors. Success in these efforts 
would significantly widen the opportunties for application of these sensors. CARBOCAP® GMT-221/222 
CO 2 sensors were analyzed in this phase of the investigation. The test team emplyed several sensors which 
were ranged from 0% to 1%, 0% to 3%, and 0% to 5%. Pre-mixed gas samples were introduced to an 
evacuated bell jar over pre-defined pressure increments. The test stand involved a bell jar, vacuum pump, 
and a series of mass flow controllers connected to pre-mixed dry gas sources, along with temperature and 
pressure sensors for data collection. 

Figure 1 demonstrates typical results obtained in this analysis for a 0-1% CARBOCAP® GMT-222 
sensor. The left panel shows the raw data collected in the investigation. The dashed line indicates the actual 
composition within the bell jar of known volume which originally starts at a given pressure and temperature. 
As the pressure in the jar is increased using the pre-mixed gas, the actual concentration within the jar 
asymptotically approaches the concentration of the gas source. As indicated in the left panel of Fig. 1, the 
reported readings deviate significantly from the actual values at low pressures. As pressure increases to 1 
atmosphere standard (atm), the compensated values and the reported readings begin to converge. 

At the time of this investigation, it became obvious compensation was necessary; however, the exact com- 
pensation relationship was not known. As a result, three potential compensation functions were considered: 
(1) compensation employing a scalar ctk to determine the degree of non-linearity associated with pressure 
compensation, (2) pressure compensation via a power-law, and (3) pressure/composition compensation via 
a power low using an equivalent pressure described by Eq. (4). The right panel of Fig. 1 demonstrates 
the results considering the three compensation functions. First, for all k measurements (on the order of 
10,000), unique ak values were determined demonstrating how significantly the reported reading needed 
re-scaled to collapse upon the actual concentration. All fitting was performed using a objective function 
and a Nelder-Mead downhill simplex method. 10 Once these values were determined, the natural logarithm 
of ak was plotted against the logarithm of pressure in atmospheres. The near linear relationship between 
the logarithms of ak and pressure seem to suggest a power law relationship is a reasonable approximation. 
Secondly, the power law compensation functions were defined along with the compensation parameters a , 
/?, </> co 2 - The values of the fitting parameters were once again determined using the simplex method and the 
relationships were plotted in both panels of Fig. 1. 

The results of this exercise are demonstrated in two particularly salient features in Fig. 1. First, for 
accurate compensation, it is clear that knowledge of both pressure and composition are necessary. Second, 
the power law relationship between pressure and reported reading means that as the pressure is decreased, 
measurement noise is dramatically amplified through the compensation algorithm. The latter observation 
suggests that COTS-NDIR sensors will have a practical limit in reduced pressure detection where the noise 
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Figure 1: Compensation of reduced pressure Vaisala CARBOCAP® GMT-221 readings for an evacuated 
bell jar re-pressurized with 0.87% C0 2 with balance of N 2 . The left panel shows sensor reported values, 
sensor values compensated according to pressure, and sensor values compensated according broadening effects 
captured by the equivalent pressure. The right panel demonstrates a unique scalar value a*,, required to 
compensate each individual reading along with how pressure and pressure/composition compensation relate 
to these readings. 


introduced through compensation renders the sensors inaccurate for application. The Vaisala sensors revealed 
that this appears to be somewhere around 0.6-0. 7 atm although it depends largely on the desired application 
of the sensor and what level of error is acceptable for test purposes. The manufacturer guarantees pressure 
compensation results in an error of 1.5% of the full-scale range plus 2% of the reading as long as the system 
pressure is between 70-130 kilopascal (kPa) (0.7-1. 3 atm) for this particular sensor. 

Results collected within this investigation demonstrated two important relationships. Sensor reported 
readings deviated significantly from actual concentration at reduced pressure and gas-broadening effects are 
significant. Since spacecraft and spacesuits are designed to provide humidity in the range of 25-75%, 11 the 
latter observation indicates that accurate compensation routines cannot be developed without accounting for 
broadening effects introduced by water vapor. Moreover, pressure relevant to spacesuit application are as low 
as 4.1 psia (0.28 atm). 6 Consequently, the next phase of the investigation sought to characterize and develop 
compensation routines for sensors capable of detection of C0 2 and H 2 0 at lower pressures. Therefore the 
Li-Cor 840A® gas analyzer was employed for this analysis that advertises a pressure compensation range of 
0.15-1.13 atm. However, during operation, the built-in pressure compensation was only found to be accurate 
when measurements deviated by approximately 0.1 atm from calibration pressure. The manufacturer later 
suggested that the supplied compensation routine should only be used at reduced pressures after re-zeroing 
and re-spanning the instrument at the pressure of interest. Since optimal pressures for spacecraft and space- 
suits is still a matter of debate, and since EVA pre-breathe protocol requires dynamic pressure acclimation 
with constant C0 2 detection, it was necessary to develop in-house compensation techniques that span the 
entire range of pressure anticipated for space exploration. 

Controlled introduction of C0 2 and water vapor was facilitated by an existing test stand described in 
detail elsewhere. 8,9 The test stand included a moisture and carbon dioxide removal swing-bed which was 
actuated to maintain an environment safe for space exploration as defined elsewhere. 11 A cavity ring-down 
mass spectroscopy (CRDS) unit provided what was considered the actual readings as discussed in a preceding 
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section. Figure 2 demonstrates the results collected during this investigation. Performance was assessed at 
four individual pressures as demonstrated by the top graph. Data were collected at the outlet side of the 
swing-bed for several cycles and then switched to the inlet of the swing- bed for several more cycles at each 
pressure. The 840A gas analyzer recorded the detection cell pressure, temperature, and reported C0 2 and 
H 2 O concentrations, thereby providing all necessary data for the compensation routine. 
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Figure 2: Compensation of reduced pressure Li-Cor 840A® readings within the swing-bed performance 
characterization test loop as compared to the Picarro CRDS unit over a variety CCU/EbOconcentrations 
and pressures. The top panel shows sensor reported values and sensor values compensated according broad- 
ening effects over multiple pressures. The bottom panel shows sensor reported values and sensor values 
compensated according broadening effects at the scale of individual cycles at 6.0 psia. 

As indicated in the top graph of Fig. 2, the reported reading differs substantially from readings collected 
by the CRDS unit before compensation. Moreover, the error appears more prominent at lower pressures 
which agrees with previous observations with other NDIR sensors. Utilizing the data collected by the 840A 
and CRDS gas analyzers, the parameters a, /3, (j ) c 02 , and </> H2 0) of Eq. (5) were fit for the 840A sensor 
via the simplex method employed in previous phases of the investigation. As evidenced in the top graph, 
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the fitted parameters resulted in good agreement with the CRDS unit. The bottom graph illustrates the 
compensated results at 6.0 psia (0.41 atm) with increased resolution at the cyclic-scale. At the resolution 
of cycles, agreement appears reasonable between the CRDS and 840A. The CRDS results demonstrate a 
‘stair-stepped’ trend which is associated with the fact that the CRDS reports to the data acquisition system 
at a frequency of 1 Hz while the measurements are actually obtained over the course of 1.6-2. 2 seconds. 
Furthermore, at the peak of each profile, an anomalous decrease in CO 2 concentration is observed. This 
artifact was introduced by the test loop which requires pneumatic actuation of the swing-bed valve when 
cycling to a fresh bed. During the actuation process, a small amount of gas was introduced into the test 
loop resulting in a short perturbation in system pressure. Whereas the CRDS has a means to control its 
detection cell pressure, the 840A does not have this capability. Consequently, in the compensation procedure, 
the pressure perturbation results in an amplification of this perturbation. In spite of these artifacts, the 
agreement between the NDIR and CRDS appears suitable. The compensated NDIR sensor results were used 
to perform a material balance on the test loop demonstrating the swing-bed removal rate and the CO 2 and 
H 2 O injection rates were in agreement. 

In summary, NDIR is an attractive method for CO 2 quantification since the technology is affordable and 
the instrumentation involved is flexible with quick installation and connection to data acquisition systems. 
However, these results in aggregate demonstrate that pressure and molecular broadening can drastically alter 
the reported readings of COTS-NDIR sensors when used for reduced pressure CO 2 detection. As demon- 
strated, these effects can be quantified and corrected. Through corrections, we were able to demonstrate 
that accurate data could be collected. This is promising for test applicability; however, the development of 
these sensors were subject to requirements that were very different than constraints imposed during space 
exploration. The latter phase of this investigation sought to develop a first principles competency in NDIR 
physics in the hopes that the knowledge could be applied in aerospace sensor design. 


V. Molecular Spectra Calculations 

Since NDIR sensors are inexpensive and flexible, it is desirable to develop methods to numerically com- 
pensate these sensors for off-nominal application. The preceding results indicate this is possible for some 
commercially available NDIR sensors which is encouraging for the terrestrial qualification of developing 
hardware. However, for application in spacecraft and spacesuits, it is not always practical to employ COTS- 
NDIR sensors. First, the sensors were not necessarily designed for these types of applications and, therefore, 
might result in noise amplification as conditions dynamically change within a closed volume. Second, the 
communication protocols from instrumentation to the core flight software of the caution and warning system 
might not be optimized. Third, the packaging of the instrumentation and processing equipment was devel- 
oped without consideration for deployment in a spacecraft and spacesuit, thereby creating problems during 
interfacing commercial counterparts with existing systems. Therefore, it is attractive to have a first prin- 
ciples understanding of how pressure, temperature, and molecular species concentrations influence infrared 
spectra. The following section outlines the prediction of spectral signatures that arise as a result of kinetic 
theory. 

Two approaches exist for predicting infrared (IR) spectra: (1) ab initio molecular dynamics calculations 
can be employed to determine the energy states and resulting IR spectra of an ensemble of molecules, 12-14 
and (2) spectral line parameters can be obtained from ‘line list’ databases such as HITRAN 15 to build spectra 
line-by-line. The former approach requires the molecular ensemble to be large for accurate predictions which 
involves a high computational cost. The latter approach can provide accuracy with limited computational 
resources. The detailed derivation for this analysis is documented thoroughly elsewhere 5 and will be explained 
briefly below. 


A. Beer-Lambert Law 


As previously discussed, absorbance at wavelength v through pathlength L is quantified according to the 
Beer-Lambert law, A = A 0 [ 1 — exp (— k v L)\. This relationship depends on the absorption coefficient, k u . 
The absorption coefficient can be written as follows for a sample that behaves as an ideal gas. 


k v 


qP 

k B T a ‘ 


( 9 ) 
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In Eq. (9), q is the volume mixing ratio (unitless), k B is the Boltzmann constant (1.380658x 10 -2 ' 5 J/K), 
and oy is the molecular cross-section (cm 2 ). The product of the component spectral line intensity S 
(cm - 1 /(molecule/cm 2 )) and spectral line shape at wavenumber u, g v ( 1 /cm -1 ), result in the molecular cross- 
section, ay = Sg v . Furthermore, noting transmissivity at v is the complement of absorption r„ = 1 — A/ 
(unitless), the following more convenient expression for the Beer-Lambert law is formulated. 


r„ = exp 


-qPL 

k B T 



(10) 


Equation (10) indicates that the absorption coefficient can be calculated for each wavelength pending the 
spectral line intensity and the line shape are known. 


B. Line Shape 

Line shape g v , is subject to a number of broadening mechanisms and is particularly sensitive to mixture 
composition and pressure. This section presumes the spectrum of a single gas in the mixture is of interest 
so that subscripts and superscripts can be neglected in the development of terms contributing to line shape 
and intensity. Subsequent sections will discuss how the species spectra are combined for the mixture trans- 
mittance spectrum. The line centers v c (cm -1 ), tend to increase linearly with pressure from a zero-center 
position v° (cm -1 ), according to the air-broadened pressure-shift coefficient 6, viz. 


V c 


+ 6 


P 

P~o 


( 11 ) 


Parameters v° and S are readily obtained for common atmospheric gases from the HITRAN database. 

While line centers is influenced linearly by pressure as captured by the pressure-shift coefficient, Doppler 
broadening arising from thermal motion blurs the line shape in a Gaussian manner providing a half-width 
Id- 


Id = v c 


2 k B T 


n l/2 


(12) 


In expression (12), rr ij is the molecular mass of component j (kg), and c is the speed of light (cm/s). In 
near-zero pressures, Doppler broadening dominates the line shape term. 

An additional species-dependent pressure-broadening effect is captured by the Lorentz line shape which 
becomes more prominent as pressure increases. In particular, collisions between similar and dissimilar 
particles result in different energy transitions. The Lorentz line width can be described by the following 
relationship. 


7 l = [(1 — q) l°L a +<H °L S . 




(13) 


In equation (13), 7 / is the air-broadened half-width contribution, is the self-broadened half- width 
contribution, and n is the coefficient of temperature dependence of air-broadened half-width. These three 
parameters can be obtained in the HITRAN database. 

The Doppler and Lorentz broadening contributions are combined yielding the Voigt profile, f v . 


fv 


Tl. f + °° ex P (~ l 2 ) 

71 - 3/2 7 f ) J ( 7i / 7 D f + [( 1 / - v c ) hn - tf 


(14) 


The function representing the Voigt profile is non-analytic. Thus a numerical technique is warranted to 
generate an approximate answer. This can be achieved through Gaussian quadrature although the infinite 
bounds must be handled carefully ( e.g . a transformation of variables can be used to modify the boundaries). 
However, the form of expression (14) lends itself to much more straightforward integration through Gauss- 
Hermite quadrature due to the numerator term which was employed in this investigation. After calculating 
the Voigt profile, one subsequent adjustment is made for far-wing effects. 


v tanh (hcv/2k B T) 
u c tanh ( hcv c /2k B T ) v 


(15) 
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In Eq. (15), h is Planck’s constant (6.62618x 10~ 34 m 2 -kg/s). It should be noted that caveats exist for special 
cases in which line-coupling and speed-dependent broadening produce important contributions. However, 
for most typical spectral calculations, the Voigt profile is considered a reasonable approximation and was 
therefore employed in this analysis. 

C. Line Intensity 

The line intensity, S, is the remaining calculation required to generate spectra for gas mixtures. The line 
intensity at off-nominal ( i.e . off-calibration) conditions requires accounting for the thermal influence on the 
rotational and vibrational partition functions, Q r and Q v . 

S(T P) = <? 6Xp ( hcEo/ksT ) 

’ ° Qr (T) Q v (T) exp {hcEo/ksTo) 

In the expression above, E 0 is the lower state energy provided by the HITRAN database. Expression (16) 
can be solved from a tabulation of total internal partition sums parameterized as a function of temperature. 16 

From this point, equations (10) through (16) provide the expressions required to generate a spectrum 
pending line and intensity parameters are known. The line and intensity parameters have been meticulously 
documented within the HITRAN database at 296 K providing all the data required to predict transmissivity 
spectra for components typically found within the atmosphere of Earth. Figure 3 illustrates the results 
of this process. In the calculations, wavenumber was discretized in increments of 0.25 cm -1 , which was 
parametrically found to be a suitable value as further discretization significantly increased computational time 
while it did not result in a different value for the integrated area under the absorptivity curve {A' = 1 — t v ). 


1 — exp (—hcv c /kBT) 
1 — exp (—hePc/ksTo) 


(16) 



Figure 3: Molecular spectrum for C02at 370 parts per million (ppm) in air. The left panel demonstrates 
the intensity of a transition in energy states over a range of wavenumber. The right panel demonstrates the 
transmissivity spectrum obtained based on intensity, broadening, and energy state data obtained from the 
HITRAN database. 15 

Figure 3(a) illustrates the line positions and intensities for carbon dioxide at 296 K. This chart shows 
strong spectral activity in the wavenumber range of 2300-2400 cm -1 . As a consequent, the Vaisala CARBOCAP® 
line is designed to detect in this range using a light source with a wavelength of 4.26 /im. The Li-Cor 840A 
provides a broad-band IR source with filters in front to provide light at 4.26 and 3.95 fan. 
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The expected IR spectrum arising from the Voigt profile calculations for 370 ppm of CO 2 in an atmosphere 
of 1 atm at 296 K is demonstrated in Fig. 3(b). The concentration of 370 ppm was chosen since it falls 
within the range of observed atmospheric CO 2 concentrations and it was the design point of the Li-Cor 840A 
sensor. An optical pathlength of 14 cm was used for the calculation which is associated with the path length 
of the 840A gas analyzer. It is evident in Fig. 3(b) that at a wavelength of 4.26 /im, the transmissivity is 
still greater than 95% for this low concentration but that it is still less than 100% indicating detection would 
be possible for a path length of 14 cm. This analysis suggests sensor design is constrained by expected CO 2 
concentration in the applied environment. Specifically, for aerospace application, the pressure, temperature, 
and expected CO 2 concentration range could be set and a parametric investigation could be performed on 
L to determine the optimal path length for the desired accuracy and resolution. It is possible that over a 
wide range of conditions, a sensor array may be required along with logic to determine which sensors from 
the array should be referenced under the measurement conditions. 


D. Transmittance Spectrum of the Gas Mixture 

Preceding sections discuss in detail the expected line shape and intensity contributions for an individual gas 
j over all spectral lines denoted by i. During an IR spectroscopy measurement of gas mixture, the observed 
profile is a coalescence of individual contributions. The prediction of the observed transmittance profile is 
generated from a product of the contributing gas spectra. 




nn^ =exp 

3 i 


PL 

k^T 


Y % Y a ^ 

3 i 


(17) 


Eq. (17) was employed for a variety of gases at different compositions in air at 1 atm and 296 K in an 
effort to verify these calculations. These results are displayed in figure 4. 



WAVELENGTH (nm) 



Figure 4: Comparison of experimentally determined spectra (left) and numerical predictions (right) for 
a mixed gas sample. Experimental results were obtained from a Vaisala product brochure for the 
CARBOCAP® sensor series. 17 


As indicated in Fig. 4, the location of IR active regions for the components agree well between the spectral 
predictions and experimentally collected results from Vaisala. 1 ' Transmissivity predictions show comparable 
agreement as well, although a path length of 14 cm was assumed since data with respect to the experimental 
instrument were not available. This analysis suggests the spectral calculations produce reasonably accurate 
results and may serve as an analysis tool for sensor design for aerospace application. Consequently, it was 
of interest to confirm source of the power law relationship discussed in previous sections through a more 
fundamental analysis. These results are demonstrated in figure 5. 

Furthermore, to explore the relationship between transmissivity and total system pressure, carbon dioxide 
concentration was fixed at 370 ppm while the system pressure was varied from 0.05 atm to 1.00 atm in 
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Figure 5: Numerical results demonstrating a loss of absorptivity as the partial pressure of air is decreased with 
constant CO 2 concentration (370 ppm), (a) The integrated area under the absorption spectrum normalized 
by the area at 1.0 atm and 296 K. (b) log-log plot of area and pressure relationship with a black-line 
demonstrating a linear fit which corresponds to compensation without accounting for composition. 


increments of 0.05 atm. Rather than compare the transmissivity predictions at 4.26 /J,m , the integrated area 
under the absorptivity curve was instead used as a surrogate in this calculation. This was chosen since it 
provides the pressure dependency across the entire IR active region of the spectrum of C0 2 . The calculated 
area was small and was subsequently normalized by the area at 1.00 atm to better demonstrate the trend. As 
shown by Fig. 5(a), the area varies non-linearly with pressure and ostensibly follows a power law relationship. 
The log-log plot in Fig. 5(b) shows that there is slight deviation from the power law relationship (solid 
black line) with pressure which supports the experimental observations displayed in Fig. 1 that component 
broadening effects are important for accurate reduced pressure C0 2 detection. Consequently, the spectral 
calculation routine seems to reproduce experimental observations and could plausibly be employed to tailor 
sensors for aerospace applications. 


VI. Conclusions 

Desired spacecraft and spacesuit pressure are still under consideration. As a result, efforts toward hard- 
ware development generally require performance characterization at a variety of pressures. COTS-NDIR 
sensors are attractive due to their flexibility and affordability for C0 2 detection. However, results demon- 
strated herein show a precipitous loss of accuracy at even modest deviations from off-calibration conditions 
which provides the impetus to develop a better understanding on NDIR theory. 

As indicated in this work, heuristic scaling arguments were utilized to develop reasonable compensation 
techniques. Experimental results confirmed this approach and provided evidence composition broadening 
significantly alters spectra when pressure is reduced. Consequently, a recursive compensation technique was 
developed with the Newton-Raphson method which was subsequently verified through experimentation. 

In an effort to better understand sensor compensation, as well as to lend insight into the design of custom 
C0 2 NDIR sensors, spectral line shape calculations were performed. The spectral calculations were verified 
with experimental results collected from literature. Moreover, the results confirm a loss of accuracy in 
reduced pressure application. The analysis seems to also provide support for the power law correction and 
the importance of accounting for composition broadening during compensation. Future work may include 
using the spectral calculation code for the parametric analysis across a variety of compositions, temperatures, 
and pressures, in the development of reduced pressure NDIR sensors or sensor arrays. 
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Mathematical Annotation 


A. Variables 

A = Absorbance, [unitless] 

A' = Relative absorbance, [unitless] 

E 0 = Lowest energy state, [J] 

F = Objective function, [unitless] 

F' = Derivative of objective function, [unitless] 

/„ = Voigt profile at u, [ 1 /cm -1 ] 

g v = Spectral luine shape at z/, [ 1 /cm -1 ] 
k v = Absorption coefficient, [m -1 ] 

L = Optical length path, [cm] 

m = Molecular mass, [kg/kmole] 

n = Iteration number, [No.] 

P = Pressure, [atm -or- psia] 

q = Volume mixing ratio, [unitless] 

Q r = Rotational partition function, [unitless] 

Q v = Vibrational partition function, [unitless] 

S = Spectral line intensity, [cm _1 /(molecule/cm 2 )] 

T = Temperature, [K] 

y = Mole fraction, [unitless] 

B. Greek Letters 

a = Compensation parameter, [unitless] 

/? = Compensation parameter, [unitless] 

S = Air-broadened pressure coefficient, [unitless] 

= Doppler-broadened half- width, [cm -1 ] 

7 l = Lorentz-broadened half- width, [cm -1 ] 

7i a = Air-broadened contribution to Lorentz half- width, [cm -1 ] 

7i s = Self-broadened contribution to Lorentz half- width, [cm -1 ] 

v = Wavenumber, [cm -1 ] 

v c = Wavenumber line center, [cm -1 ] 

v° = Wavenumber line center at reference conditions, [cm -1 ] 

(j) = Broadening coefficient, [unitless] 

<j v = Molecular cross-section, [m 2 ] 

r = Convergence tolerance, [unitless] 

r„ = Transmissivity at v, [unitless] 
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C. Sub/super-scripts 

0 = Associated with reference conditions 
v = Associated with a given wavenumber 

1 = Denotes a particular species 
j = Denotes a particular species 

D. Physical Constants 

c = Speed of light, [2.997925x10 s m/s] 
h = Planck’s constant, [6.626180x 10 -34 J-s] 
ks = Boltzmanns constant, [1.380658x 10 -23 J/K] 
P 0 = Reference pressure, [1 atm -or- 14.696 psia]] 
R = Universal gas constant, [8.314472 J/mol-K] 
T 0 = Reference pressure, [296.15 K] 
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